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We establish an intriguing connection between quantum phase transitions and bifurcations in the reduced 
fidelity between two different reduced density matrices for quantum lattice many-body systems with symmetry- 
breaking orders. Our finding is based on the observation that, in the conventional Landau-Ginzburg- Wilson 
paradigm, a quantum system undergoing a phase transition is characterized in terms of spontaneous symmetry 
breaking that is captured by a local order parameter, which in turn results in an essential change of the reduced 
density matrix in the symmetry-broken phase. Two quantum systems on an infinite lattice in one spatial dimen- 
sion, i.e., quantum Ising model in a transverse magnetic field and quantum spin 1/2 XYX model in an external 
magnetic field, are considered in the context of the tensor network algorithm based on the matrix product state 
representation. 
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Introduction. Recently we have witnessed a growing in- 
terest in the study of quantum many-body systems in the 
context of the fidelity approach to quantum phase transi- 
tions (QPTs) |1, 2] with_symmetry-breaking/topological or- 
ders \XA.5.§.Z&.^^Mm- InRefs. 14^43, it 
has been argued that the ground state fidelity per site may be 
used to detect QPTs. Since the argument is solely based on 
the basic Postulate of Quantum Mechanics on quantum mea- 
surements, it is expected that this approach is applicable to 
quantum lattice systems in any spatial dimensions, regardless 
of what type of internal order is present in quantum many- 
body states. In fact, it has been confirmed that the ground 
state fideUtyper site is able to describe QPTs arising from an 
SSB lIlElMSH, the Kosteriitz-Thouless transition [II and 
topological QPTs in the Kitaev model [10]. Remarkably, the 
ground state fidelity per site may be computed in terms of the 
newly-developed tensor network (TN) a lgorithms, such as the 
matrix product states (MPS) ||T3, 14, 15] in one spatial dimen- 
sion, and the tensor product states (TPS) Olg] , or equivalently, 
the projected entangled-pair states (PEPS) jlTll . in two and 
higher spatial dimensions. 

In the conventional Landau-Ginzburg-Wilson paradigm, an 
SSB, which occurs when a system possesses a certain symme- 
try whereas the ground state wave functions do not preserve 
1911 . is quantified in terms of a local order parameter. 
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An intriguing feature of local order parameters for quantum 
systems with symmetry-breaking orders is bifurcations they 
exhibit at critical points. Remarkably, such a bifurcation also 
manifests itself in the ground state fidelity per site [8]. The ad- 
vantage of the latter over local order parameters lies in the fact 
that the ground state fidelity per site is universal in the sense 
that it is not model-dependent, in contrast to model-dependent 
order parameters in characterizing QPTs in quantum lattice 
many-body systems. 

The investigation that has been caiTied out so far mainly 
focuses on the ground state wave functions for quantum lat- 
tice systems. Thus the fidelity is equivalent to the overlap 
between two different ground states, which are pure states. 



However, as first discussed in Ref. fd\ (see also [fl2l|), the 
fidelity defined for two different mixed states is also useful 
to locate phase transition points. The mixed states are de- 
scribed by the reduced density matrices arising from tracing 
out the degrees of freedom in the environment surrounding a 
sub-system. Indeed, for a sub-system of a composite quan- 
tum system, a reduced density matrix is a basic notion that 
is indispensable for the analysis of a composite quantum sys- 
tem. In fact, the presence of a local order parameter makes 
a difference in the reduced density matrices in the symmet- 
ric and symmetry-broken phases. From this one may antici- 
pate that the (reduced) fidelity between two different reduced 
density matrices for quantum lattice many-body systems with 
symmetry-breaking orders should capture bifurcations arising 
from an SSB. 

In this paper, we attempt to address this problem. We shall 
investigate two quantum models on an infinite lattice in one 
spatial dimension, i.e., quantum Ising model in a transverse 
magnetic field and spin 1/2 XYX model in an external mag- 
netic field. Both systems possess a discrete symmetry group 
Z2, but the ground states break the symmetry. The SSB is 
reflected as a bifurcation ["s*] in the ground state reduced fi- 
delity for both systems. Our result demonstrate that one may 
identify a phase transition point as a bifurcation point in the 
ground state reduced fidelity between two different reduced 
density matrices for quantum lattice many-body systems with 
symmetry -breaking orders [20]. 

The models. The first model we consider in this paper is 
the quantum Ising model in a transverse magnetic field in an 
infinite-size lattice in one spatial dimension. The Hamiltonian 
takes the form. 



(1) 



where S a {a - x,z) are the spin 1/2 Pauli operators at site /, 
and A is the transverse magnetic field. The model is invariant 
under the symmetry operation: ^i'' — > and 5''^ 5''^ 



2 



for all sites, which yields the Z2 symmetry. As is well known, 
the system undergoes a second order QPT at the critical field 

= 1 im. 

The second model is the spin 1/2 XYX model in an external 
magnetic field. The Hamiltonian can be written as 

H^Yj + Ay5;,'l5|,'+'l + 5 W^fi J + hS^\ (2) 

/'=— 00 

where S J^' {a - x, y, z) are the Pauli spin operators at site /, A,, 
is a parameter describing the anisotropy in the internal space, 
and h is the external magnetic field. This model also possesses 
a Z2 symmetry, with the symmetry operation: 5^'' S^!^, 
-S\^ and 5''' 5''^ for all sites. Below we shall 
choose Ay = 0.25. In this case, the critical magnetic field is 
h, ~ 3.210(6) I2I. 

The reduced density matrix. We are now in a position to 
clarify the difference of the reduced density matrices in the 
symmetric and symmetry-broken phases for both quantum 
Ising model in a transverse magnetic field and quantum spin 
1/2 XYX model in an external magnetic field. The analysis 
will be carried out for both the one-site and two-site reduced 
density matrices, respectively. For the quantum Ising model 
in a transverse magnetic field, the one-site reduced density 
matrix in the Z2 symmetric phase takes the form, 

Pising^^+2{S,)S„ (3) 

where (S^) is the expectation value of 5 ^ in the ground state in 
the Z2 symmetric phase, whereas the two-site reduced density 
matrix in the Z2 symmetric phase is. 

Pi sing = - / + 4q' 1 5 X ® 5 X + 4ff 2^5 ^ ® 5 J 

+ 03! IS> S , + rlSl I. (4) 

Here, ai = {Sx^S_^}, 02 - {S. ®S-), a^, - (I^S.), and 
a4 - {S^^ I), with / being the identity matrix. 

In the Z2 symmetry-broken phase, the presence of the 
nonzero local order parameter {S x) implies that the one-site 
reduced density matrix takes the form, 

Pising^^+2{S.,)Sx + 2{S,}S,. (5) 

Such a violation of the symmetry is also reflected in the two- 
site reduced density matrix: 

Pising - — I + ^a\S jc ® S x + 4-028 c ® 5 j 

+4a6Sz®Sx + aTl<8iSx + asSx<S>I, (6) 

with 05 = (5x ® Sz), a(, - {S , ® 5.t), ai - {I ® S x), and 
0-8 = (5 ,-«)/). 

For the quantum spin 1/2 XYX model in an external mag- 
netic field, the one-site reduced density matrix in the Z2 sym- 
metric phase takes the form, 

PXYX^\+2{S,}S„ (7) 



while the two-site reduced density matrix is 

PxYx = i/4-4/?i5,®5, +4/?25y®5y 

+4/3-iS:(E)S,+/34l®S,+/35S,,'S)I, (8) 

with ySi = {Sx ® Sx), P2 = (Sy ® 5v), yS3 = (S, ® 5 J, /34 = 
</®5,), and/?5 = (5;®/). 

In the symmetry-broken phase, the one-site reduced density 
matrix becomes 

PxYX^\+2{Sx}Sx + 2{S,}S„ (9) 

whereas the two-site reduced density matrix is 

pxYx = + 4^iSx®Sx+4/32Sy(»Sy +4/338-^8, 

+l}4l®8z+ 1358^®! + 4l3(,8x®8z (10) 
+4l3i8, ®8x+l3iI®8x+li98x®I, 

with Pf, = {8x ® 8,), Pi = {8, ® 8x), ySg = (/ ® 8 x), and 
ySg = {8x®I). 

Bifurcations in the reduced fidelity between two reduced 
density matrices and quantum phase transitions. The reduced 
fidelity measures the distance between two quantum mixed 
states. Specifically, for two reduced density matrices and 
Pi', the reduced fidelity F{p_\,px) is defined to be 

F(pA,Px)^tr^pfp,.pf. (11) 

Here, p 1 and px are the reduced density matrices correspond- 
ing to two diff'erent values, A and X , of the control parameter 
A. Notice that the reduced fidelity F(pA,px) is a function of A 
and A , which satisfies the following properties: (i) normaliza- 
tion F(pa,Pa) = 1; (ii) symmetry F(pA,p/) = F(px^Pa)\ (iii) 
range 0< F(pa,Pa') <1- 

The connection between a bifurcation point in the ground 
state partial fidelity between two reduced density matrices and 
a critical point for quantum lattice systems undergoing QTPs 
with symmetry-breaking order may be established in the fol- 
lowing way. In the conventional Landau-Ginzburg-Wilson 
paradigm, a quantum system undergoing a QPT is charac- 
terized in terms of an SSB that is captured by a local order 
parameter, which in turn results in an essential change of the 
reduced density matrix in the symmetry broken phase, as seen 
above. More precisely, consider a quantum system on an in- 
finite lattice in one spacial dimension, with Z2 as a symme- 
try group [23]. Suppose the Z2 symmetry is spontaneously 
broken when the control parameter A crosses a critical point 
A = Ac, then the system undergoes a QPT with a nonzero local 
order parameter in the symmetry-broken phase. Such a local 
order parameter is defined in a local area Q on a lattice, and 
it is not invariant under the action of the symmetry operation 
generating the symmetry group Z2, thus yielding degenerate 
ground states in the symmetry-broken phase. Therefore, the 
reduced fidelity F(pA,px), with the reference state p^ in the 
symmetry-broken phase, yields two different values, which 
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correspond to two distinct values of the local order parame- 
ter arising from two degenerate ground states, if is in the 
symmetry-broken phase, and one value, if is in the sym- 
metric phase. That is, a bifurcation point occurs in the ground 
state reduced fidelity between two reduced density matrices 
F(pa,P/), with the critical point being the bifurcation point. 
However, if we choose the reference state p ,' in the symmetric 
phase, then no such a bifurcation occurs, due to the invariance 
of the reduced density matrix p/ under the action of the sym- 
metry operation. In addition, the same argument applies to 
the reduced fidelity defined on any local area Q' with f2 as a 
subset [24]. 

The results. In the context of the TN algorithm initi- 
ated by Vidal [25], the problem to find the system's ground 
state wave functions amounts to computing the imaginary 
time evolution for a given initial state |*I'(0)): \^{t)) - 
exp(-//T)|>I'(0)>/|exp(-//T)|>I'(0))|. An efficient way to 
achieve this task is to exploit the Suzuki-Trotter decomposi- 
tion [26], which allows us to reduce the imaginary time evo- 
lution operation to a product of two-site evolution operators 
acting on sites / and / -i- 1: U{i,i + 1) = exp(-/z''''^'-'^i5T), 
6r « 1. In addition, any wave function admits an MPS 
representation in a canonical form: attached to each site is 
a three-index tensor or and to each bond a diago- 
nal singular value matrix Aa or Ag, depending on the evenness 
and oddness of the i-th site and the i-th bond. Here, 5 is a 
physical index, s = I, - ■ ■ ,d, with d being the dimension of 
the local Hilbert space, and I and r denote the bond indices, 
l,r- 1, - ■ ■ ,x, with X being the truncation dimension. The 
action of a two-site gate U{i, i + 1) may be absorbed by per- 
forming a singular value decomposition, thus resulting in the 
update of the MPS representation. Repeating this procedure 
until the ground state energy converges, one may generate the 
ground state wave functions in the MPS representation. We 
emphasize that, in practice, we adjust the truncation dimen- 
sion;^ to identify a critical point from bifurcation points in the 
reduced fidelity F{p,i,px)- In fact, it saves a lot of the compu- 
tational resources to perform simulations for relatively small 
values of x- Usually, a shift in the bifurcation points occurs 
due to the finiteness ofx- Therefore, it is necessary to perform 
an extrapolation with respect io x to locate the critical point. 

In Figdia), we plot the ground state reduced fidelity 
F{p,\,px) between the one-site reduced density matrices for 
the quantum Ising model in a transverse field with the field 
strength A as the control parameter Here, we choose px {X - 
0.9) as a reference state, which breaks the Z2 symmetry. 
The one-site reduced fidelity can distinguish two degenerate 
ground states with a bifurcation point as a pseudo phase tran- 
sition point A^ [8?]. When the parameter A is tuned beyond 
such a pseudo transition point, two degenerate ground states 
vanish, implying that the system undergoes a phase transition. 
One observes that, with increasing, the pseudo phase transi- 
tion point A^ moves toward the exact value 1 . Performing an 
extrapolation of A^ with respect to;^', we get Ac - 1.00233. In 
FigHfb), we show the two-site reduced fidelity for the quan- 
tum Ising model in a transverse magnetic field. The same ref- 
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FIG. 1 : (color online) (a) The ground state one-site reduced fidelity 
F(Pa,Pa' ) for the quantum Ising model in a transverse magnetic field. 
We choose p^' (with A = 0.9) as the reference state, which is in 
the Z2 symmetry-broken phase. The pseudo phase transition point 
/l^ occurs as a bifurcation point. When we enlarge the truncation 
dimension x, the pseudo phase transition point /i^ is getting closer 
to the exact value. Inset: the critical point A,, is determined from 
an extrapolation of the pseudo phase transition point with respect 
to the truncation dimension x- Here, the fitting function is A^ = 
A, + a^-*, with Ac = 1.00233, a = 0.39373 and b = 2.43904. (b) The 
ground state two-site reduced fidelity, F(pi,p^'), for the quantum 
Ising model in a transverse magnetic field. The same reference state 
p/ (A = 0.9) has been chosen as in the case of the one-site reduced 
fidelity. Both insets indicate the same pseudo critical points for the 
ground state one-site and two-site reduced density matrices for the 
system. This is expected due to the fact that they are resulted from 
the same set of the ground states. 



erence state is selected as in the case of the one-site reduced 
fidelity. We observe that a bifurcation also occurs in the two- 
site partial fidelity. Indeed, it yields the same pseudo phase 
transition points A^, and thus the same critical point Ac- 

In Fig. |2a), the ground state one-site reduced fidelity 
F{ph,Ph' ) for the quantum spin 1/2 XYX model in an external 
magnetic field h is plotted, with the magnetic field strength h 
as the control parameter We choose py, with hi - 3.05, in 
the Z2 symmetry-broken phase, as the reference state. The re- 
duced fidelity F(pi,,pi/) distinguishes two degenerate ground 
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FIG. 2: (color online) (a) The ground state one-site reduced fidelity 
F{ph,Ph' ) for the quantum spin 1/2 XYX model in an external mag- 
netic field. Here, the external magnetic field strength h is the control 
parameter. We choose p,,' ( with h = 3.05) as the reference state, 
which is in the symmetry-broken phase. The pseudo phase transi- 
tion point occurs as a bifurcation point. For a larger value of 
the pseudo phase transition point is getting closer to the known 
value 3.2049. Inset: the critical point he is determined from an ex- 
trapolation of the pseudo phase transition point with respect to the 
truncation dimension The fitting function is = he + a^"*, with 
K = 3.2049, a = 0.0267 and b = 0.9426. (b) The ground state two- 
site reduced fidelity F{pt„p^' ) for the quantum spin 1/2 XYX model. 
We still choose p,/ {h = 3.05) as the reference state. The two-site 
reduced fidelity exhibits the similar behavior to that of the one-site 
reduced fidelity. Both insets display the same pseudo critical points 
for the one-site and two-site reduced fidelity for the system. 



states in the symmetry-broken phase. The bifurcation point 
resuhed from the one-site reduced fideHty is the pseudo 
phase transition point, which is quite close to the known crit- 
ical value he ~ 3.210(6) f22"], if the truncation dimension^ is 
large enough. Performing an extrapolation with respect to x 
yields the critical point he - 3.2049. In Fig. |2jb), we plot 
the two-site reduced fidelity for the quantum XYX model. 
Here, the same reference state as in the case of the one-site 
partial fidelity has been chosen. The two-site reduced fidelity 
F(ph,ph') is also able to detect the Z2 SSB. Indeed, it again 
yields the same bifurcation points /i^, and thus the same criti- 



cal point he- 

Summary. We have established an intriguing connection 
between QPTs and bifurcations in the reduced fidelity be- 
tween two different reduced density matrices for quantum lat- 
tice many-body systems with symmetry-breaking orders. Our 
work is based on the newly-developed TN algorithms, which 
produce degenerate ground states in the symmetry-broken 
phase for quantum lattice systems under QPTs arising from 
an SSB. Two quantum systems on an infinite lattice in one 
spatial dimension, i.e., quantum Ising model in a transverse 
magnetic field and quantum spin 1/2 XYX model in an exter- 
nal magnetic field, have been investigated in the context of the 
TN algorithm based on the MPS representation. 
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